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Dependence of transition frequencies on the fine structure constant a = e 2 /he is calculated for 
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I. INTRODUCTION 

Theories unifying gravity with other interactions as 
well as many cosmological models allow for space-time 
variation of fundamental constants. Experimental search 
for the manifestation of this variation spans the whole 
lifetime of the Universe from Big Bang nucleosynthesis 
to the present-day veryprecise atomic clock experiments 
(see, e.g. reviews [l|, Q). An evidence that the fine- 
structure constant might be smaller about ten billion 
years ago was found in the analysis of quasar absorp- 
tion spectra @, H H H 0, @] • This finding together with 
progress in developing of very precise atomic frequency 
standards motivated many laboratory searches for the 
present-day time variation of the fundamental constants 
(see, e.g. [9(). In particular, strong limit on the rate 
of the time variation of the fine structure constant a 
(a = e 2 /he) were found by comparing frequencies of dif- 
ferent atomic transitions over few years Q. 

Apart from the microwave atomic clocks and optical 
frequency standards, a number of atomic transitions in 
which the change of frequency due to change of a is 
strongly enhanced has been suggested in Refs. [hJ HH, 

Interpretation and planning of the measurements of 
the a variation require atomic calculations to relate the 
change of atomic frequencies to the change of the fine 
structure constant. A number of such calculations for 
atomic optical transitions have been performed in our 
early works P3, EH, 03, 03, 03, [E|. Independent cal- 
culations for some optical transitions have been recently 
reported in Ref. 

From the computational point of view the most impor- 
tant parameter of an atom which determines the choice 
of the computational method as well as the accuracy 
which can be achieved in the calculations is the num- 
ber of electrons in open shells. The larger the number 
the more difficult are the calculations. Many optical fre- 
quency standards are based on atoms or ions with just 
one or two valence electrons (9|. Calculations for such 
systems are accurate and reliable p3, 03 03] • How- 
ever, many atomic systems which are used or planned 
to be used in laboratory search for variation of the fine 



structure constant have more then ten electrons in open 
shells. For example, strong limits on the variation of al- 
pha in time[l7j and variation of alpha due to change of 
the gravitation potential fl8| were obtained with the use 
of dysprosium atom which has twelve external electrons 
(see also Q3, 03, 0). There are plans to use holmium 
(13 electrons) for similar measurements [20|. There are 
ongoing measurements or plans for measurements for 
Ag I !2r], Yb II, Jp H, Yb III 24] and Hg II [H 
(see also a review Q and references therein). These sys- 
tems involve states with excitations from d or / subshells 
and therefore must be treated as many-valence-electrons 
systems. 

Calculations for many-valence-electron atoms are dif- 
ficult due to the fast growth of the matrix size of the 
configuration interaction (CI) eigenvalue problem with 
the increase of the single-electron basis. In our recent 
paper on Fe l[2(| we used a version of the CI method 
which is similar to multi-configuration CI method (see, 
e.g. Ref. [13]) and which allows to obtain reasonably ac- 
curate result with a very short basis. In present paper 
we use this method for many-electron systems which are 
of the interest for laboratory search of the variation of 
the fine structure constant. The aim of the calculation is 
to check our early results as well as to calculate relativis- 
tic energy shifts for atomic transitions which have never 
been considered before. 



II. METHOD 

Detailed discussion of the method can be found in our 
early works [ToL [26j. Here we repeat its major points. 

It is convenient to present the dependence of atomic 
frequencies on the fine-structure constant a in the vicin- 
ity of its physical value ao in the form 

u(x) =u> + Qx, (1) 

where ujq is the present laboratory value of the frequency 
and x = (a/ao) 2 — 1, g is the coefficient which is to be 
found from atomic calculations. Note that 

9 = -j- • (2) 

dx x=0 
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To calculate this derivative numerically we use 

^ uj(+S) - u)(-S) 
? ~ 2S ' 



(3) 



In the present calculations we use S = 0.05, which leads 
to 



10 (w(+0.05) - w(-0.05)) . 



(4) 



In a single-electron approximation relativistic energy 
shift can be estimated using the formula [l(| 



1 



3a + 1/2 



— C{Z,j a , l a ) 



(5) 



where a is the index for a single-electron state, E a is 
its energy, v a is its effective principal quantum number 
(y a = \/\/—2E a ), j a and l a are total and angular mo- 
menta of the state a. C(Z,j a ,l a ) is a parameter which 
is introduced to simulate the effect of Hartree-Fock ex- 
change interaction and other many-body effects. For a 
transition between many-electron states which can be ap- 
proximated as a single-electron transition from state a in 
lower level to state b in upper level one has 



q 1 



(6) 



The formulas ©,([1]) are too inaccurate for practical use 
in the interpretation of the measurements. However, they 
are very useful for predicting what one can expect to 
find in different atomic transitions and for explaining the 
values and sign of the relativistic corrections. We will 
use it for the discussion of our results. 

For accurate numerical calculations of the coefficients q 
using ((3]), a must be varied in the computer code. There- 
fore, it is convenient to use a form of the single electron 
wave function in which the dependence on a is explicitly 
shown (we use atomic units in which e = Ti = 1, a = 1/c) 



njlm 



1 / f v (r)Q(n)ji m 
r \ iag v (r)il(n)ji ri 



(7) 



where n is the principal quantum number and an index 
v replaces the three- number set n,j,l. This leads to a 
form of radial equation for single-electron orbitals which 
also explicitly depends on a: 



df v 
dr 



+ —fv(r)- 2 + a 2 (e v -V HF ) g v (r) = 



dg v K v . . . 

-j fv{r) + [e v 

dr r 



V H F)fv(r) = 0, 



(8) 



here K = (-l) l+j+1/2 (j + 1/2), and V HF is th e Har tree- 
Fock potential. Equation ([5]) with a — a^^f 5 + 1 and 
different Hartree-Fock potential Vhf for different config- 
urations is used to construct single-electron orbitals. 

Table U lists configurations considered in present work. 
For Ag I, Au I and Hg II we use only ground-state config- 
uration and configurations, involving excitation from the 
upper core <i-state. The latter corresponds to the states 



TABLE I: Configurations and effective core polarizabilities 
(a p , a.u.) used in the calculations. 



Atom 


Z 


N v a 


Set 


Parity 


Configuration 


a p 


Agl 


47 


11 


1 


Even 


4d lu 5s 


0.4 








2 


Even 


4d 9 5s 2 


0.414 


Dy I 


66 


12 


1 


Even 


4/ 10 6s 2 


0.4 








2 


Even 


4/ 10 5d6s 


0.397 








3 


Even 


4/ 9 6s 2 6p 


0.4039 








4 


Even 


4/ 9 5d6s6p 


0.389 








5 


Odd 


4/ 9 5d 2 6s 


0.3895 








6 


Odd 


4/ 9 5d6s 2 


0.4 








7 


Odd 


4/ 10 6s6p 


0.393 


Ho I 


67 


13 


1 


Odd 


4/ n 6 S 2 


0.4 








2 


Odd 


4/ 10 6s 2 6p 


0.401 








3 


Odd 


4/ n 5d6s 


0.401 








4 


Odd 


4/ 10 5d6s6p 


0.39 








5 


Odd 


4/ n 6p 2 


0.39 








6 


Even 


4/ 10 5d6s 2 


0.3927 








7 


Even 


4/ n 6s6p 


0.3962 








8 


Even 


4/ 10 5d 2 6s 


0.39 








9 


Even 


4/ 10 5d6p 2 


0.4 








10 


Even 


4/ 10 6s6p 2 


0.4 


Yb II 


70 


15 


1 


Even 


4/ 14 6s 


0.4 








2 


Even 


4/ 13 6s 2 


0.399 








3 


Even 


4/ 13 5d6s 


0.3911 








4 


Even 


4/ 13 5rf 2 


0.39 


Yb III 


70 


14 


1 


Even 


4/ 14 


0.4 








2 


Even 


4/ 13 5d 


0.3914 








3 


Even 


4/ 13 6s 


0.3977 


Au I 


79 


11 


1 


Even 


5d 10 6s 


0.4 








2 


Even 


5d 9 6s 2 


0.417 


Hg II 


80 


11 


1 


Even 


5d 10 6s 


0.4 








2 


Even 


5d 9 6s 2 


0.426 



a N v is the number of valence electrons. 



which are to be used in the measurements. We add more 
configurations for Yb II and Yb III and even more for 
Dy I and Ho I. In the latter atoms the states of interest 
are highly excited ones for which configuration mixing is 
strong and should be taken into account more accurately. 

The self-consistent Hartree-Fock procedure is done for 
every configuration listed in Table U separately. Then va- 
lence states found in the Hartree-Fock calculations are 
used as basis states for the CI calculations. It is impor- 
tant for the CI method that the atomic core remains the 
same for all configurations. We use the core which cor- 
responds to the ground state configuration. Change in 
the core due to change of the valence state is small and 
can be neglected. This is because core states are not sen- 
sitive to the potential from the electrons which are on 
large distances (like 6s, 6p and hd electrons). The 4/ 
electrons are on smaller distances and have larger effect 
on atomic core. However, in all the cases (see Table HJ 
only one among about ten 4/ electrons change its state. 
Therefore their effect on atomic core is also small. More 
detailed discussion of the effect of valence electrons on 
atomic core can be found in Refs. 128. 12911. 
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The effective Hamiltonian for N v valence electrons has 
the form 

N v n v 
H^ = ^2h u + Y^e 2 / rij , (9) 

i—l i<j 

hi(n) is the one-electron part of the Hamiltonian 

Ze 2 

hi = ca ■ p + (8 - l)mc 2 + V core + SV. (10) 

r 

Here a and 8 are Dirac matrixes, V core is Hartree-Fock 
potential due to core electrons and SV is the term which 
simulates the effect of the correlations between core and 
valence electrons. It is often called polarization potential 
and has the form 



Here a p is polarization of the core and a is a cut-off 
parameter (we use a = ag). 

The form of the SV is chosen to coincide with the stan- 
dard polarization potential on large distances (— a v /2r ). 
However we use it on distances where valence electrons 
are localized. This distances are not large, especially 
for the 4/ electrons. Therefore we consider SV as only 
rough approximation to real correlation interaction be- 
tween core and valence electrons and treat a p as fitting 
parameters. The values of a p for each configuration of 
interest are presented in Table HI They are chosen to fit 
the experimental position of the configurations relative 
to each other. For all configurations of the same atom 
the values of a p are very close. This is not a surprise 
since the core is nearly the same for all configurations of 
interest. One can probably say that small difference in 
a p for different configurations simulates the effect of in- 
completeness of the basis and other imperfections in the 
calculations. 

Tables [H] and IIIII present comparison between exper- 
imental and theoretical energies and <?-factors for Dy I 
and Ho I atoms. The ^-factors are useful for the iden- 
tification of the states and for control of configuration 
mixing [30] . For dysprosium atom both the energies and 
g-factors are reproduced quite accurately. This includes 
the states with the energies of 19797.96 cm -1 which are 
used in the measurements (T3, EH • 

For holmium the g-factors are not known for the most 
of the states. This makes it more difficult to identify the 
states and to judge about the accuracy of the calculations 
of the relativistic energy shifts. If the measurements for 
holmium are to go ahead it would be good to measure 
the g-factors as well, at least for the states of most in- 
terest. At the moment we can only rely on the energies. 
Although the energies are reproduced in the calculations 
quite accurately the coefficients q in ([1]) are very sensitive 
to the configuration mixing which in turn is sensitive to 
the energy intervals between close levels of the same par- 
ity and total angular momentum. Therefore, having good 



TABLE II: Energy levels (cm ) and g-factors of some low 
states of Dy I 



Conf. 


J 


Experiment 


a 


Calculations 










9 




9 


4/ iu 6s 2 


8 


u.uu 


i 
i 


.242 


u 


1 O/IOQ 




7 


A 1 1A 9*} 


i 

i 


.173 


A A HQ 


1 1 7A7 

1.1/4/ 




6 


i UOU.Ol 


i 
i 


.072 


( ouu 


l.U / Zo 




5 


Q91 1 PiS 
i?Z 1 l.Oo 


n 
u 


.911 


iJiJOO 


u.yuou 




4 


luyzo.zo 


n 
u 


.618 


1 1 a/in 

llo4U 


U.DlDo 


4/ 10 5d6s 


9 


i 7p;i a 

1 ( 014. OU 


i 

i 


.316 


1 / / Uo 


1.O140 




8 


loyUo.zl 


1 


.22 


lyooD 






7 


91 (Y7 A 9fl 
Z1U / 4.ZU 


i 


.24 


91 QQ1 


l.lyoo 




8 


1 / Olo.OO 


i 


.33 


1 7S71 
1(0(1 


i oofin 




7 




I 


.28 




1.3012 




6 


91 1 7Q 


I 


.24 


22042 


1.2116 




7 


18094 V) 


I 


.38 


18308 


1.3835 




6 








20090 


1.3078 




5 








22478 


1.2198 




10 


18462.65 


1 


.282 


18461 


1.2883 




9 


19240.82 


1 


.217 


19592 


1.2277 




8 


20193.60 


1 


.16 


20893 


1.1700 


4/ 10 6s 2 


8 


19019.15 


1 


.14 


21377 


1.1113 


4/ 10 5d6s 


11 


19348.72 


1 


.27 


19295 


1.2675 




10 


1 Q7Q7 QR b 


1 


.21 


20077 


1.2089 




9 


20209.00 


1 


.14 


20847 


1.1261 


4/ 9 6s 2 6p 


7 


20614.32 


1 


.32 


19835 


1.3372 




8 


20789.85 


1 


.32 


19832 


1.2997 


4/ 10 5d6s 


8 


21603.04 


1 


.26 


23205 


1.2514 




7 


21778.43 


1 


.26 


23232 


1.2419 




9 


22045.79 


1 


.22 


23429 


1.2677 




10 


22487.14 


1 


.197 


24132 


1.2162 


4/ 9 5d6s6p 


8 


23031.46 


1 


.37 


23132 


1.3730 


4/ 9 5d6s 2 


10 


12892.76 


1 


.29 


12920 


1.2933 


4/ 10 6s6p 


10 


17513.33 


1 


.30 


17582 


1.2944 


4/ 9 5d 2 6s 


10 


19797.96'' 


1 


.367 


19693 


1.3677 


4/ 9 5d 2 6 


10 


21788.93 


1 


.34 


22312 


1.3340 



"NIST, Ref. [3ij 

''States used in the measurements [l7, 18] 



accuracy for absolute values of energies is not enough for 
reliable results for the coefficients q. It is very important 
that the relative positions of the states around the states 
of interest are reproduced accurately in the calculations. 



III. RESULTS AND DISCUSSION 

a. Holmium Holmium atom has been suggested for 
the search of the variation of the fine structure constant 
by Mark Saffman [20j. From the computational point of 
view it represents the most difficult case. It has thirteen 
electrons in open shells, very dense spectrum, strong con- 
figuration mixing and multiple level crossing in the vicin- 
ity of the physical value of a when energies are considered 
as functions of a 2 . All these factors contribute to the in- 
stability of the results. Therefore, it is instructive to 
start from simple estimations based on a single-electron 
approximation. Table IIVI shows approximate values of 
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TABLE III: Energy levels (cm 1 ) and g- factors of some low 
states of Ho I 



Conf. 


Parity 


J 


Expt." 


Calculations 








Energy 


Energy 


g 


4f n 6s 2 


Odd 


15/2 


0.00 





1.20 






13/2 


5419.70 


5770 


1.11 


4/ 10 6s 2 6p 


Odd 


15/2 


18572.28 


18343 


1.28 


4/ n 5d6s 


Odd 


13/2 


18867.40 


18684 


1.37 






15/2 


19276.94 


19295 


1.32 


4/ 10 5d6s6p 


Odd 


15/2 


24112.04 


23908 


1.35 


4/ 10 5d6s 2 


Even 


15/2 


8427.11 


8395 


1.28 






13/2 


9147.08 


9341 


1.33 


4/ 10 5d6s 2 


Even 


15/2 


12339.04 


12903 


1.23 






13/2 


12344.55 


12953 


1.23 






11/2 


13082.93 


13799 


1.25 






13/2 


15081.12 


16459 


1.17 






11/2 


16937.43 


16817 


1.13 


4/ n 6s6p 


Even 


15/2 


15855.28 


15913 


1.28 






13/2 


17059.35 


17135 


1.20 


4/ 10 5d 2 6s 


Even 


15/2 


20167.17 


20138 


1.41 



NIST, Ref. [H 



TABLE IV: Approximate values of the (/-coefficients for dif- 
ferent configurations of Ho I (xlO 3 cm -1 ). 



Configuration 


Parity 


Transition" 


q 


4/ n 6s 2 


Odd 


(ground state) 


5(5) 


4/ 10 5d6s 2 


Even 


4/^5d 


-35(15) 


4/ n 6s6p 


Even 


6s — > 6p 


4(4) 


4/ 10 6s 2 6p 


Odd 


4/^6p 


-45(15) 


4/ n 5d6s 


Odd 


6s — > 5d 


7(4) 



"single-electron transition from the ground state. 



q-coefficicnts for different configurations of holmium ob- 
tained with the use of formula (O but with energy shifts 
of the individual single-electron states (A a ^) taken from 
the Hartree-Fock calculations rather than from formula 
©. Note that we present the energies and ^-coefficients 
relative to the ground state. Therefore, relativistic en- 
ergy shift (q) for the ground state is zero by definition. 
The ^-coefficients for other states of the 4/ n 6s 2 config- 
uration are determined by the fine structure of the 4/ 
orbital. Large error bars are due to the fact that rela- 
tivistic energy shifts depend on the values of the total 
momentum j of the single-electron states involved in the 
transition (see formula ©). For example, the 6s — > 6p 
transition can be the 6s — > 6pi/2 or the 6s — > 6P3/2 tran- 
sition, etc. 

As can be seen from the data in Table IfVI the values 
of q are very different for different configurations. For 
example, q w —35000cm -1 for the 4f w 5d6s 2 configura- 
tion and q w 4000cm -1 for the 4/ 11 6s6p configuration. 
But these two configurations have the same parity and 
can have the states of the same total angular momentum 
J. Therefore, if these configurations are strongly mixed 



the resulting values of q will be linear combinations of 
q « —35000cm -1 and q w 4000cm -1 , i.e. they may take 
any value between large negative value and some positive 
one depending on which configuration dominates in the 
state. The same is true for odd states which constitute a 
mixture of the negatively shifted 4/ 10 6s 2 6p configuration 
with the positively shifted 4/ 11 5d6s or 4/ 11 6s 2 configu- 
rations. The analysis of the holmium spectrum shows 
that there are many states of the same parity and to- 
tal angular momentum J which are separated by only 
small energy intervals and in which different configura- 
tions dominate. These states are strongly mixed which 
leads to instability of the calculations of the g-coefficients. 
The only way to obtain reliable results is to make sure 
that the relative position of the states in the vicinity of 
state of interest as well as the energy intervals between 
these states are reproduced accurately in the calculations. 
This can be achieved by appropriate choice of the a p pa- 
rameters for different configurations (see Table UJ). 

In Table fVl we present the results of the calculations for 
two pairs of almost degenerate states of holmium. The 
relative change of frequency between degenerate levels 
due to change of a is strongly enhanced by small en- 
ergy interval. The enhancement factor K (defined by 
Slu/uj = KSa/a where u> is the transition frequency) is 
given by [l(| 

K = 2Aq/AE. (12) 

This enhancement factor is about 3 x 10 5 for both pairs 
of holmium states presented in Table fVl To avoid misun- 
derstanding we should note that the enhancement of the 
relative effect here is due to the small AE; there is no 
any enhancement of the absolute values of the frequency 
shifts. The values of q in holmium are typical for heavy 
atoms. 

b. Dysprosium Dysprosium atom is used for the 
search of time-variation of the fine structure constant at 
Berkeley [H EI El. It has two almost degenerate lev- 
els of opposite parity at E = 19797.96cm -1 for which 
the enhancement factor (fT2"]) is about 10 s [T3|. Limits 
on the rate of changing of alpha in time obtained from 
monitoring the frequency of the transition between these 
two levels over long period of time is on the same level of 
precision as for the most advances atomic optical clock 
experiments (~ 10 -15 /yr) [l7j]. The interpretation of the 
measurements are based on our early calculations [l4j]. 
The aim of present calculations is to check our previous 
result with a significantly different method. 

From the computational point of view dysprosium 
atom is an easier case than holmium in two ways. First, it 
has one less valence electron and its spectrum is much less 
dense. The energy separation between mixing states is 
larger than 1000 cm -1 which is much easier to reproduce 
in the calculations than few hundred cm -1 as in the case 
of holmium. Second, experimental values for g-factors 
are available for dysprosium. The ^-factors are almost as 
sensitive to configuration mixing as the q-coefficients [30l ] 
providing an important test of the accuracy of the calcu- 
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TABLE V: Experimental and theoretical energies and calculated relativistic energy shifts ((/-coefficients, cm 1 ) for some tran- 
sitions of Ag I, Dy I, Ho I, Yb II, Yb III, Au I and Hg II. 



Atom 


Ground state 


Excited state 


Energy [cm ] 


g-coefficients, 


[cm" 1 ] 




Conf. 


J 


Conf. 


J 


Expt." 


this work 


this work 


other 


Agl 


4d iu 5 S 


1/2 


4d y 5s 2 


5/2 


30242.26 


30188 


-11300 










4d 9 5s 2 


3/2 


34714.16 


35114 


-6500 




Dyl 


4/ 10 6s 2 


8 


4/ 10 5d6s 


10 


19797.96 


20077 


7952 


6008 6 








4/ 9 5d 2 6s 


10 


19797.96 


19693 


-25216 


-23708 6 


Ho I 


4/ 11 6s 2 


15/2 


4/ 10 5d6s 2 


11/2 


20493.40 


21763 


-28200 










4/ 11 5d6s 


13/2 


20493.77 


20872 


7300 










4/ n 6s6p 


13/2 


22157.86 


22599 


3000 










4/ n 5d6s 


9/2 


22157.88 


22631 


8000 




Yb II 


4/ 14 6s 


1/2 


4/ 13 6s 2 


7/2 


21418.75 


20060 


-63752 


-56737'' 








4/ 13 6s 2 


5/2 


31568.08 


31303 


-53400 










4/ 13 5d6s 


5/2 


26759.02 


26781 


-46863 










4/ 13 6s6p 


7/2 


47921.31 


47927 


-60432 




Yb III 







4/ 13 5d 





45276.85 


46505 


-32800 


-27800 6 


Au I 


5d 10 6s 


1/2 


5d 9 6s 2 


5/2 


9161.3 


9186 


-38550 










5d 9 6s 2 


3/2 


21435.3 


22224 


-26760 




Hg II 


5d 10 6s 


1/2 


5d 9 6s 2 


5/2 


35514.624 


35066 


-52200 


-56670 c 








5d 9 6s 2 


3/2 


50555.567 


50886 


-37700 


-44000 c 



a NIST, Ref. |31| 
6 Dzuba et al, Ref. 
c Dzuba et al Ref. 



lation. As can be seen from Table |TT] both energies and 
g-factors are reproduced in the calculations with good 
accuracy. 

The results for the ^-coefficients are compared in Ta- 
ble [V] with previous calculations. The largest relative 
difference is for the smaller coefficient and is about 30%. 
However, the difference for Aq which is important for the 
interpretation of the measurements is only 12%. 

c. Ag I, Yb II, Yb III, Au I and Hg II The Ag I, 
Yb II, Yb III and Hg II atoms are also used or considered 
for the use in the laboratory search for the variation of the 
fine structure constant (see @ and references therein). 
We have included Au I because it has electron structure 
similar to Ag I and Hg II. However, we are unaware about 
any plans to use Au in the measurements. 

All these systems utilize the use of a transition from 
the ground state to a low lying state which involves an 
excitation from the core. Both states have no significant 
admixture of other configurations, relatively easy to cal- 
culate and produce stable results. 

The results for the g-cocfficicnts for the transitions of 
interest are presented in TablefV] Here we also have good 



agreement with previous calculations for the cases when 
the data are available. 



IV. CONCLUSION 

Calculations of the relativistic energy shifts are pre- 
sented for many transitions in many-valencc-electrons 
systems which are used or planned to be used in the labo- 
ratory search for variation of the fine structure constant. 
Good agreement with previous calculations confirms the 
analysis based on old results and provides an estimate 
of the accuracy of the calculations. Many atomic transi- 
tions are added which were never considered before. 
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